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Abstract 

It is argued that the dominant contribution to the interaction of 
quark gluon plasma at moderate T > T c is given by the nonperturba- 
tive vacuum field correlators. Basing on that nonperturbative equa- 
tion of state of quark-gluon plasma is computed and in the lowest 
approximation expressed in terms of absolute values of Polyakov lines 
for quarks and gluons Lf un d(T); L a dj(T) = (Lj un d) 9 / 4 known from lat- 
tice and analytic calculations. Phase transition at any [i is described 
as a transition due to vanishing of one of correlators, D (x), which 
implies the change of gluonic condensate AG2. Resulting transition 
temperature T c (/i) is calculated in terms of AG2 and Lf un( i{T c ). 

The phase curve T c (fi) is in good agreement with lattice data. 
In particular T c (0) = 0.27; 0.19; 0.17 GeV for n f = 0,2,3 and fixed 
AG 2 = 0.0035 GeV 4 . 



1 Introduction 

QCD at nonzero temperature T and quark chemical potential /1, and specifi- 
cally phase transitions are important topics both for theory, and experiment 
[H[2] and are carefully studied on the lattice (for recent reviews see [31 E]). 

Since these phenomena cannot be explained in perturbative QCD, any 
analytic approach should be based on nonperturbative methods, and as such 
solvable models, like Nambu-Jona-Lasinio model, have been used mostly. 
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However the models used do not contain confinement and treat phase transi- 
tion from the chiral symmetry (CS) point of view. Therefore it is interesting 
to look at the Equation of State (EoS) and phase transition in the framework 
of a nonperturbative (NP) method, based on vacuum fields and implementing 
confinement. 

In this letter we are using such tool, which is called the Field Correla- 
tor Method (FCM) In FCM all NP information about QCD vacuum is 
contained in field correlators, and the Casimir scaling, now proved on the 
lattice also for T > T c [6], allows ( with ~ 1% accuracy) to consider only 
four quadratic correlators for T, \i > of colorelectric (colormagnetic) field 
strengths, D E (x), Df (x), (D H (x), D±(x)). Here D E ,Df are defined as vac- 
uum averages of colorelectric field Ej, 



and similarly for D ,D± with colormagnetic fields Hi instead of E i} where 
the sign of the last term is reversed, and y) = Pexpig J A^dz^ is par- 
allel transporter (Schwinger line). Note that only D E contributes to the 
string tension a E implying confinement, and D H to the spatial string ten- 
sion, o H = a s , while Df, contain both perturbative and nonperturbative 
parts. 

The correlators were calculated on the lattice [H [8] and also analytically 



Most calculations with FCM at T — \i = contain as input only string 
tension a, strong coupling a s (q) and current (pole) quark masses, and results 
for hadronic spectra are in a good agreement with available experimental and 
lattice data, see last ref. in [5] and [10] for reviews. 

In [11] the FCM was extended to the case of nonzero T, and in [12] 
the phase transition at \i — was described as a Vacuum Dominated (VD) 
transition from the confining; vacuum (all D^,Df\i = E,H are nonzero, 
and D E = D H , Df = Df 1 as for T = 0) to the deconfined state (only D E 
and hence a E = \ f D E (x)d 2 x are zero). This picture was supported by 
lattice measurements [3 [8] where indeed D E was shown to vanish at T = T c , 
while D E and colormagnetic correlators stay intact. This fact supports the 
conjecture [13] that D E is nonzero at T > T c and strong enough to support 
bound states of quarks and gluons. Bound states of qq were indeed found on 



g 2 (tr f [Ei{x)${x,y)E k {y)<$>{y,x)}) 

QD E ~\ , 
5 ik D E + D E + ul—± +u iUk - 




dDf 



du 2 



[9j. 
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the lattice years later (see [3] for a review) and in [TH [To] the connection of 
those to the correlator Df was quantitatively established. 

Moreover in [H] it was shown, that a new static potential Vi(r, T) is 
provided by Df at T > T c and the NP part of it yields the modulus of the 

renormalized Polyakov line, Lf un d = exp (— Vl ^.' T ' > j ■ 

The physical picture of the VD deconfining phase transition given in [12j , 
has allowed to get the (simplified) estimate of the transition temperature T c 
at n — in reasonable agreement with lattice data for rif = 0, 2, 3 the only 
input being the standard value of the gluonic condensate G2, known from the 
ITEP sum rule method [16] . It is the purpose of the present letter to extend 
the method to the case of nonzero fi and nonzero interaction with vacuum 
(Vi 7^ 0) and to find T c (/z), for /i > 0. 

To this end we are using the EoS found recently in [17] , where pressure in 
the lowest approximation is expressed through the only dynamical quantity 
- LfundiT) calculated via the correlator Df (x) or obtained from lattice data. 
This factor takes into account interaction of the quark with the NP vacuum in 
the form of Lf un d(T) and analogous factor for gluons is L a dj(T) = (Lf un d) 9 ^. 
The interaction between quarks and gluons is considered as a next step, 
and is argued to possibly contribute in a narrow temperature region near T c 
[T7| IT8] . This simple picture of VD dynamics with the only input Lf und (T), 
taken from lattice or analytic calculation suggested in [17] . which may be 
called Vacuum Dominance Model (VDM), is adopted below and is shown 
to produce surprisingly reasonable results, being in good agreement with 
available lattice data for /i = or nonzero /1, where these data are reliable. 



2 Nonperturbative EoS for /i > 

The main idea of the VDM discussed below is that the most important part 
of quark and gluon dynamics in the strongly interacting plasma is the inter- 
action of each individual quark or gluon with vacuum fields. This interaction 
is derived from field correlators and is rigorously proved to be embodied in 
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factors, which happen to coincide with the modulus of Polyakov loopS, 

£/w = exp( — ), L adi = exp(-- — ) (1) 

where Vi(T) = Vi(oo,T), V D = V D (r*,T) and V^T),^ found in [H] to 
be (/? = 1/T) 

Vi(r,T) = f du{l-uT) f ^Df(^eT^), (2) 
Jo Jo 

Vb(r,T) = 2 f dv{l-vT) f\r - ^D E (y/e + ^). (3) 
Jo Jo 

In Vd, Eq.([T]) r* is the average size of the heavy- light Qq or its adjoint 
equivalent system; for T > T c one has D E = Vd = 0. For T < T c one has 
r /und = oo for n/ = 0, yielding Lf un d = 0, however r* dj - m 0.4 fm for any 
n/ and gives nonzero L a dj(T < T c ). The form ([TJ)- (J3j) is in good agreement 
with lattice data [6] and also explains why Lf und is a good order parameter 
(however approximate for rif ^ 0). Note, that only NP parts D E ,D E enter 
in d2J), dSD, see first ref. in [9] for discussion of separation of these parts; 
renormalization procedure is discussed in [TH [15] . 

In the lowest NP approximation one neglects pair, triple etc. interactions 
between quarks and gluons (which are important for T c < T < 1.2T C where 
density is low and screening by medium is not yet operating, see pTTj [18])) 
and derives the following EoS (this approximation is called in [17] the Single 
Line Approximation (SLA)) 

Pg ~ T 4 ~ n 2 n 4 ^fundVq C0Sn j, l 4 J 



n=l 



p,f A _ 2(iv c 2 - r 



Vgl 



T 4 7T 2 ^ n 

n=l 



v 

^4 L ^i ( 5 ) 



with 



^(T)= n A(^f)^-\(^y + ... (6, 



: We neglect in the approximation the difference between Lf unc i expressed via Vi(oo, T) 
and L^nd where the role of V\(T) is played [HI [17] by singlet QQ free energy F^q(oo, T). 
The latter quantity contains all excited states, so that Vi(T) > Fi~(oo, T). 
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<C -L — 



-|x|/A 



1 GeV, where A is the vacuum 
hence powers of Lf, see [T7j for 



In (J3J, ©it was assumed that T 
correlation length, e.g. d[ e \x) ~ e" 
details. 

With few percent accuracy one can replace the sum in (J5j) by the first 
term, n — 1, and this form will be used below for while for p q this re- 
placement is not valid for large ^, and one can use instead the form equivalent 
to (SD, 



Pi 



7T 2 



where v = m q /T and 



T 



(7) 



1 



7^ 



3 V 2 2 +i/ 2 — C 



(8) 



Eqs. ([7]), ([5]) define p q ,p g i for all T, /i and m q , which is the current (pole) 
quark mass at the scale of the order of T. 

To draw p q ,p g i and p = p q + p g i as functions of T, fi one needs explicit 
form of V\(T). This was obtained analytically and discussed in [14] ; another 
form was found from Df (x) measured on the lattice in [HI [9] and is given in 

m. 

Also from lattice correlator studies [HE], P3] Vi(T = T c ) is (with ~ 10% 
accuracy) 0.5 GeV and is decreasing with the growth of T (cf. Fig. 2 of [14] 
and Fig. 1 of [IS]). This behaviour is similar to that found repeatedly on the 
lattice direct measurement of see e.g. Fig. 2 of [19] where = V X {T) 
is given for nj = 0,2,3. In what follows we shall exploit the latter curves 
parametrizing them for T >T C and all nj as 



Vi(T) 



0.175 GeV 



1.35 ( jr 



V X {T C ) w 0.5 GeV. 



(9) 



For n > one can expect a /i-dependence of Vi, however it should be weak 
for values of fi much smaller than the scale of change of vacuum fields. The 
latter scale can be identified with the dilaton mass ma, which is of the order 
of the lowest glueball mass, i.e. ~ 1.5 GeV = m d . Hence one can expect, that 
V\ in the lowest approximation does not depend on fi. This is supported by 



the lattice measurements in 



where for 



1.5 and ^ 



0.8 the values 



of are almost indistinguishable from the case of ^ 



0. 
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Figure 1: Pressure ^ from Eq. (j4"f5]) as function of temperature T (in GeV) 
for n/ = 3, 2, (top to bottom) and AG 2 = 0.0034 GeV 4 . 

To give an illustration of the resulting EoS we draw in Fig.l the pressure 
p for the cases fi = 0, rif = 0,2. One can see a reasonable behaviour for 
T c (n f ) = (0.27; 0.19; 0.17) GeV (for n f = 0,2,3 respectively) similar to the 
lattice data, see [21 J for a review and discussion. 

3 Phase transition for nonzero fi 

Here we extend the VD mechanism suggested in [12j to the case of nonzero 
/i and V\. We assume as was said in Introduction that the phase transition 
occurs from the full confining vacuum with all correlators D E , Df, D H , 
present to the deconfined vacuum where D E vanishes. The basics of our 
physical picture is that all fields (and correlators) do not change for /i, T 
changing in a wide interval, unless /i, T become comparable to the dila- 
ton mass m& ~ M(glueball + ) ~ 1.5 GeV(E Therefore correlators and 
a E , a H are almost constant till T = T c and at T > T c a new vacuum phase 

2 Note that G*2 does not depend on fi, n./ in the leading order of the 1/N C expansion and 
one expects a growth of the magnetic part of G2 at T > Tdim.red ~ 2T C , G™ ag w 0(T A ), 
in the regime of dimensional reduction. 
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with D E = a E = is realised,which yields lower thermodynamic potential 
(higher pressure). Lattice measurements [TJ [S] support this picture. The 
crucial step is that one should take into account in the free energy F of 
the system also the free energy of the vacuum, i.e. vacuum energy density 

Evac = \0w = j^}(( F ^) 2 ) = - {11 ~ 3 l nf) G2, which can be estimated via the 
standard gluonic condensate [S] G 2 = ^((Fp 2 ) « 0.012 GeV 4 . 

Hence for the pressure P = —F one can write in the phase I (confined) 

Pi = \e vac \+Xi(T) (10) 

where x(T) is the hadronic gas pressure, starting with pions, Xpion — ^T 4 . 
In the deconfined phase one can write 

Pn = \e d v e a c c\ + (P 9 i + P q )T 4 (11) 

where \s^ c c \ is the vacuum energy density in the deconfined phase, which is 
mostly (apart from D E (0) ~ 0.2D E (0) see [7] and first ref. of [H]) colormag- 
netic energy density and by the same reasoning as before we take it as for 
T = 0, i.e. |£ mc | = 0.5|£ mc |. 

Equalizing Pj and Pu at T = T c (fi) one obtains the equation for T c 



V Pgi +Vq J 



(12) 



where A|e„ ac | = \e vac \ - |e^| w (n 3 | n/) AG 2 ; AG 2 ~ \G 2 ] p g i and p q are 
given in (jSJ), © respectively and depend on both T c and //. 

In this letter we shall consider the simplest case when the contribution 
of hadronic gas Xi(T) can be neglected in the first approximation. Indeed, 
pionic gas yields only ~ 7% correction to the numerator of (TT2"j) at T w T c , 
and from [22] one concludes that x(T c ) < 0.5T C 4 , which yields a < 10% 
increase of T c for G 2 ~ 0.01 GeV. 

From the expression for T c (/i) ([121 one can find limiting behaviour of 
T c (yU — » 0) and /i c (T — > 0). For the first one can use for p q and p g (jlj) and 
(J5j) and expand r.h.s. of ( |T2|) in ratio with the result. 
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where the last term yields a 3% correction, and = ^ n 
Solving ([TBI for T c one has 



r < - ^ m ( 1 + t 1 + i^r) (14) 

with k = |Vi(T c ). From fTl2|) . ([TBI one can compute expansion T c (yu) in 
powers of /x, 

T e Qi B ) = T c (0) (l - <7jg|y) , /is = 3//. 



C = 1 + ^ + T(0) . = 0.0110(3) for 7i* = 2, 3, 4. 

One can see, that C practically does not depend on rif and is in the same 
ballpark as the values found by lattice calculations, see [23] , |3] for reviews 
and references. 

Another end point of the phase curve, /i c (T — > 0), is found from ([T2l 
when one takes into account asymptotics $o( a ~ * °°) = \ + ^ft 2 + + •••> 
which yields (for small 



niq 

i 1 



P.AT 0) = ^ + (48) 1 "T(») (l + 



v 2 (--^)V 



(15) 

The resulting curve T c (/x) according to Eq.( fl2l with xi = is given in 
Fig. 2 forAG 2 = 0.00341GeV 4 , n/ = 2, 3 and m 9 = 0. 



4 Discussion of results 

Our prediction of T c (/i) depends only on two numbers: 1) the value of gluonic 
condensate AG 2 and 2) the value of Vi(T c ) = 0.5 GeV taken from lattice data 
[T9] (and quantitatively close to the value from the analytic form [14]). 

We take G 2 in the limits 0.004 GeV 4 < G 2 < 0.015 GeV 4 , the value 
G 2 = 0.008GeV 4 (AG 2 = 0.0034 GeV 4 ) being in agreement with lattice data 
of T c (0), for rif = 0, 2, 3, see Table 1. 



8 



0.1 0.2 0.3 0.4 0.5 0.6^ 



Figure 2: The phase transition curve T c (/x) from Eq. (jl2]) (in GeV) as function 
of quark chemical potential fi (in GeV) for rif = 2 (upper curve ) and rif = 3 
(lower curve) and AG 2 = 0.0034 GeV 4 . 



Table 1. The values of T c (/i = 0) and /x c (T = 0) computed using (TH"|) and 
( TO)]) for several values of AG 2 and = 0, 2, 3 
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AG 2 
0.01 GeV 4 


0.191 


0.341 


0.57 


1 


T c ( GeV) n f = 


0.246 


0.273 


0.298 


0.328 


T c ( GeV) n f = 2 


0.168 


0.19 


0.21 


0.236 


T c ( GeV) = 3 


0.154 


0.172 


0.191 


0.214 


/x c ( GeV) n f = 2 


0.576 


0.626 


0.68 


0.742 


/i c ( GeV) n./ = 3 


0.539 


0.581 


0.629 


0.686 



Note that T c at rif = in Table 1 is obtained not from (jUj) . but directly 
from ([HD with n/ = 0, x(T) = 0. 

The curve in Fig. 2 has the expected form, which agrees with the curve, 
obtained in |24j by the reweighting technic and agrees for [i < 300 MeV with 
that, obtained by the density of state method [23], and by the imaginary fi 
method [26]. 

An analysis of the integral ([8]) for v = reveals that it has a mild singular 
point at ix S i n g — y ± inT, which may show up in derivatives in At 
T = 0, Using = -y — 0.25 GeV and is close to the point where one expects 
irregularities on the phase diagram [23J [27] . 

The limit of small fi is given in fll4p . Taking Vi(T c ) w 0.5 GeV as follows 
from lattice and analytic estimates, one obtains with ~ 3% accuracy the 
values of T c given on Fig. 2 for nj = 2, 3 and AG2 = 0.0034 GeV 4 . 

The values of /z c , from f[T5|) are given in Table 1 and are in agreement 
with the curves shown in Fig.3. Note that x(T = 0) = and (fT31) holds also 
in the case, when hadron (and baryon) gas is taken into account. 

At this point one should stress that our calculation in VDM of p(T) does 
not contain model parameters and the only approximation is the neglect 
of interparticle interaction as compared to the interaction of each one with 
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the vacuum (apart from neglect of x{T))- Fig.l demonstrates that this ap- 
proximation is reasonably good and one expects some 10-i-15% accuracy in 
prediction of T c (/x). Note that in VDM the phase transition is of the first 
order, which is supported for n/ = by lattice data, see e.g. |6j however 
for n/ = 2, 3 lattice results disagree, see [281 [29] for a possible preference of 
the first order transition. One however should have in mind that the final 
conclusion for lattice data depends on input quark masses and continuum 
limit. 

The "weakening" of the phase transition for n/ > and nonzero quark 
masses is explained in our approach by the flattening of the curve P(T) at 
T pa T c when hadronic gas xCO is taken into account, since x(T) for rif = 
is much smaller than for rif = 2. 

A few words about chiral symmetry properties of the transition. It was 
argued in a series of papers of one of the authors, that Chiral Symmetry 
Breaking (CSB) and confinement are closely connected, and moreover the 
effective scalar quark mass operator was calculated in terms of D E (x)(a E ) 
in [SU] and chiral condensate (qq) and f n via o E in [3T] . 

Therefore (qq) and f n disappear at the same T c where D E (x) (and a E ) 
vanishes. This fact is in perfect agreement with known lattice data and 
supports the results of the present study. 

Finally, as seen from our expressions for p q , p g (jlj), (J7|), our EoS is indepen- 
dent of Z(N C ) factors and Z(N C ) symmetry is irrelevant for the NP dynamics 
in our approach. This result is a consequence of more general property - the 
gauge invariance of the partition function for all T, fi, which requires that 
only closed Wilson loops appear in the resulting expressions, yielding finally 
only absolute value of the Polyakov loop, or in other words, is expressed via 
only singlet free energy of quark and antiquark F^q ~ V±. 

5 Summary and conclusions 

We have calculated the QCD phase diagram in the lowest order of the back- 
ground perturbation theory taking into account the change of gluonic vac- 
uum condensate in the phase transition and modification of individual quark 
and gluon propagators in vacuum fields. Some support of this independent 
single particle approach can be found in the lack of strangeness (S')-charge 
(baryonic B or electric Q) correlations on the lattice [32] for T > 1.1T C . 
The resulting phase curve T c (/i) depends only on two fundamental param- 
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eters: 1) the change of gluonic condensate across phase boundary, AG 2 ~ 
and 2) the absolute value of Polyakov loop at T = T c ,Lf und (T c ) = 

exp \ — ^zT^ ) ' or e ^ se singlet Qq f ree energy V%(T C ) = F^q(oo,T). 

The second quantity is known from lattice and analytic calculations, 
V X {T C ) w 0.50(5) GeV. ' 

It is shown that T c (/i) depends on AG2 rather mildly and yields reasonable 
values of T c (0) « 0.2 GeV for generally accepted G 2 ~ 0.01 GeV 4 with 
/i c (T = 0) » 0.6 GeV. Moreover for the same value of AG 2 = 0.00341 GeV 4 
one obtains the set of three values T c (0) = (0.27; 0.19; 0.17) GeV 4 all being 
in agreement with lattice data. 

Our results do not contain model or fitting parameters and allow in 
this (however crude) approximation connect in a very simple way the phase 
boundary to the fundamental properties of the QCD vacuum. 

Many additional quantities like quark number susceptibilities, baryon - 
strangeness correlation etc. can be easily calculated in the method and shall 
be given elsewhere, as well as influence of possible diquark pairing. 

The authors are grateful for discussions to N.O.Agasian, K.G.Boreskov, 
B.L.Ioffe, A.B.Kaidalov, O.V.Kancheli, B.O.Kerbikov. 
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